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We define the entropy operator as the negative of the logarithm of the density matrix, give a 
prescription for extracting its thermodynamically measurable part, and discuss its dynamics. For 
an isolated system we derive the first, second and third laws of thermodynamics. For weakly-coupled 
subsystems of an isolated system, an expression for the long time limit of the expectation value of 
the rate of change of the thermodynamically measurable part of the entropy operator is derived 
and interpreted in terms of entropy production and entropy transport terms. The interpretation 
is justified by comparison to the known expression for the entropy production in an aged classical 
Markovian system with Gaussian fluctuations and by a calculation of the current-induced entropy 
production in a conductor with electron-phonon scattering. 


I. INTRODUCTION 

Attempts to show how nonequilibrium thermodynamic 
behavior emerges from the underlying quantum mechan¬ 
ics of individual particles is now being dubbed quan¬ 
tum thermodynamicsii^— Several approaches have arisen, 
revealing important aspects in this endeavor, such as 
how thermal fluctuations and external driving mecha¬ 
nisms affect the stochastic course of nonequilibrium pro¬ 
cesses of small systems;^ which has led to fluctuation 
theorems^^— going beyond the results from the Kubo lin¬ 
ear response theory, as well as generalized fluctuation- 
dissipation relations as studied in isolated quantum sys¬ 
tems after a quenchi^^— Other aspects, more in the spirit 
of traditional nonequilibrium statistical mechanics,— in¬ 
clude thermalization in isolated quantum systems^^— 
and the establishment of steady states in open quan¬ 
tum systems— A unified treatment along the lines 
of the classical theory of nonequilibrium thermodynam¬ 
ics is of crucial importance for a clear identification of 
the quantum-to-classical correspondence and the new 
features brought about by fully quantum-mechanical 
nonequilibrium behavior. 

The remarkable success of the classical theory2^“— in 
the description of macroscopic phenomena in fluids moti¬ 
vates us to ask what are the basic ingredients of this for¬ 
malism that such a unified treatment of quantum thermo¬ 
dynamics must also contain. We remind that the building 
blocks of the classical theory are: (i) macroscopic observ¬ 
ables, explicitly defined as a set of thermodynamically 
measurable or slowly-varying quantities, (ii) conservation 
laws for these variables and, as a foundational pillar, (iii) 
an entropy balance equation is established, splitting the 
rate of change of entropy as a part which is irreversibly 
produced, in accordance with the second law of thermo¬ 
dynamics, and a part which is transported. The validity 
of this theory relies on the local equilibrium assumption, 
whereby the nonequilibrium thermodynamic entropy is 
considered locally as a function of the same extensive 
variables as in equilibrium. 

Although significant attempts to give a meaning to 
entropy out of equilibrium's have long been known in 
quantum statisticSfSSr— a complete theory of quantum 


entropy production has not been provided yet. The main 
problem is how to conceive an adequate quantum entropy 
balance equation without assuming local equilibrium. 

For an isolated system, there is no entropy to be trans¬ 
ported outside the system and hence the entropy balance 
equation reduces to finding the right quantum expression 
for entropy whose rate of change is non-negative, accord¬ 
ing to the second law of thermodynamics, this rate then 
being the entropy production. Important efforts have 
been devoted to obtain such an expression from the den¬ 
sity matrixfiSiSSiSS but the third law of thermodynamics, 
involving the vanishing entropy of pure states has not 
been satisfactorily established. 

On the other hand, for a subsystem of an isolated 
system the establishment of a quantum entropy balance 
equation has been partially addressed^^— by assuming 
that the rate of change of an adopted expression for the 
nonequilibrium entropy of the subsystem, obtained from 
the reduced density matrix, is directly connected, as in 
the classical theory, with the rate of change of its energy. 
This involves the identification of a microscopic expres¬ 
sion for heat which is not unique^ and therefore quite 
problematic, but most importantly, does not constitute a 
full deviation from the local equilibrium assumption, as 
we show later. 

The purpose of this paper is to provide a more gen¬ 
eral treatment of quantum entropy production and then 
lay the foundation of a unified theory of quantum ther¬ 
modynamics in close correspondence with the classical 
theory. We introduce a new thermodynamic entropy op¬ 
erator St for isolated quantum many-body systems and 
show that the rate of change of its expectation value is 
non-negative, according to the second law of thermody¬ 
namics. Unlike previous approaches, we establish the 
third law of thermodynamics as a well-defined vanishing 
of the thermodynamic entropy for pure states. 

The quantum entropy balance equation for a given sub¬ 
system of an isolated system is obtained by first study¬ 
ing the time evolution of (dtSt) for the isolated system 
from first principles, i.e. from the Liouville-von Neu¬ 
mann equation for the density matrix, using the stan¬ 
dard generalized master equation approach of nonequi¬ 
librium statistical mechanics^^— and by subsequently 
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making reasonable assumptions regarding the factoriza¬ 
tion properties of the nonequilibrium probability distri¬ 
bution of microscopic states over the degrees of freedom 
of the different subsystems. 

We restrict here to weakly-coupled subsystems to show 
how our theory is consistent with the classical theory, to 
elucidate the manner in which the local equilibrium ap¬ 
proximation can be fully abandoned, and to pave the way 
to study cases of strong coupling between subsystems for 
which the aforementioned factorization properties of the 
probability distribution of microscopic states become the 
main subject of study, marking a deep connection with 
quantum information theory. A detailed investigation of 
a new methodology to approach these cases will be con¬ 
sidered elsewhere. 

The pursue of the so outlined research program is es¬ 
sential both for a more fundamental understanding of 
nonequilibrium behavior^ and because entropy produc¬ 
tion is inherent to dissipation so that a good atomic-scale 
description may have technological impact, e.g. by en¬ 
abling better control of waste heat and thermoelectric 
effects in single-molecule electronicsguiding the ef¬ 
ficient design of quantum refrigerators^ and quantum 
heat machines^ nanosized photoelectric devices^^ nan¬ 
othermoelectric engines^iiS based on quantum dots, etc., 
which are envisioned as practical applications of quantum 
thermodynamics. 

It turns out, as we show here with a particular example 
of electronic conduction in the presence of phonon modes 
playing the role of a reservoir, that our theory gives an 
explicit expression for the Joule heating from a calcu¬ 
lation of the steady state electronic entropy production 
alone. This represents an important progress since this is 
done without calculating the rate of change of the energy 
of the electron subsystem. 

This paper is organized as follows: in section |TT] we 
give a brief review of entropy production and the second 
law of thermodynamics as they manifest in phenomeno¬ 
logical thermodynamics. In section IIIII we discuss the 
local equilibrium assumption from a quantum perspec¬ 
tive, with a derivation of the first law of thermodynamics 
from the expression for (dtSt) in this case, which is shown 
to hold for quasistatic transformations or slow processes. 
This section, which mainly discusses how the foundations 
of the classical theory are to be understood quantum- 
mechanically, serves as a motivation to introduce the op¬ 
erator dtSt for general isolated quantum systems of which 
possible reservoirs are part of. 

A transition is made in section m to the generalized 
thermodynamic description of quantum systems. The 
second and third law of thermodynamics are established 
here for any isolated system, and an entropy balance 
equation is derived, splitting (dtS) into entropy pro¬ 
duction and entropy transport terms. In section El we 
show how the theory is consistent with Onsager’s clas¬ 
sical stochastic entropy production in an aged system. 
Finally, in section IVTl we calculate the electronic entropy 
production in a simple metal consisting of independent 


electrons weakly coupled to phonons in the presence of 
an external electric field, deriving the Joule heating, and 
we conclude with section IVIIl 

II. ENTROPY PRODUCTION IN 
PHENOMENOLOGICAL THERMODYNAMICS 

The thermodynamic definition of entropy changes for 
any kind of process in a closed system (not interchanging 
particles with the reservoirs) was given by Clausius at 
the very end of his monumental 1865 paperi ^^:^° If the 
system, which is considered to be in contact with a set of 
heat sources at different temperatures T, follows a path 7 
in the space of thermodynamic states, joining the initial 
and final arbitrary states A and B, respectively, then the 
thermodynamic entropy change in the process is 

Sb-Sa = Nch] + r{dQ/T)^, (1) 

J A 

where dQ is an infinitesimal amount of heat absorbed 
from (or surrendered to) the heat source at temperature 
T, and the quantity Ycjy], representing what came to 
be known as the “uncompensated heat of Clausius” ,— is 
a functional of the process. Clausius defined it in such a 
way that 

{dQ/T)^ - [ idQ/T)^,, 

JB 

, , , ( 2 ) 

where y/j is an arbitrary reversible path which is “imag¬ 
ined” to bring the system back to its initial state A. He 
proved that 

> 0, (Clausius’ inequality), (3) 

for any 7, which was a generalization of Carnot’s results 
for cyclic processes; the equality holding if and only if 
7 is a reversible path. This is the starting point of all 
the discussions found in textbooks of the second law of 
thermodynamics;^ and is therefore regarded here as the 
fundamental expression for this law. 

A classical formulation of nonequilibrium thermody¬ 
namics has been founde d^^i^^ by taking as starting point 
0 written in differential form and generalized to apply 
locally in small volume elements, Ju, of a system 

dS = diS + deS, (4) 

where diS = dNc is the entropy produced, during an in¬ 
finitesimal time interval, due to irreversible processes tak¬ 
ing place inside the volume element, and deS the entropy 
supplied from its surroundings (= dQ/T for a closed ele¬ 
ment). The second law of thermodynamics requires only 
that the entropy produced satisfies 

diS > 0, (Clausius’ inequality). (5) 

The theory so obtained for the phenomenological entropy 
production, H^^ = diS/dt, successfully describes slow 


Nell] = -i dQ/T = - 
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processes or phenomena where the decay time of local 
perturbations is very short compared to the global relax¬ 
ation time, as in chemical reactions, diffusion processes, 
heat conduction, and their cross effects in gases and liq¬ 
uids. However, it requires fundamental modifications for 
fast processes2Zi2^ and, in the following, we argue from 
a quantum-mechanical perspective why this happens to 
be the case, setting the stage and motivating the method 
for the subsequent development of our theory. 

III. LOCAL EQUILIBRIUM AND QUASISTATIC 
QUANTUM TRANSFORMATIONS 

Consider an isolated macroscopic system, possibly con¬ 
taining a set of particle and heat reservoirs which is di¬ 
vided into macroscopic subsystems. Microscopically, the 
total system is defined by the Hamiltonian 

H = E Hi -f ^ ^ Him, (6) 

l l<m 

where Hi is the Hamiltonian of subsystem I, involving 
the kinetic energies of the particles comprising the sub¬ 
system as well as the energy of interaction among all 
these particles; and Him is the Hamiltonian respresent- 
ing the interactions among the particles of subsytem I 
with those of subsystem m, possibly including hopping 
terms allowing particle transfer. 

The fundamental assumption of statistical mechanics^ 
is that, since the interaction energy among the parts 
scales with their common surface areas, while the en¬ 
ergy of the parts scales with their respective volumes, 
we can then remove all Him in from a macroscopic 
description of the dynamics and introduce instead a set 
of time-dependent parameters embodying macro¬ 

scopic constraints for the subsystem I, that evolve in 
time due to changes in the other subsystems. The op¬ 
erator representing macroscopic energy measurements in 
this approximation is 

n = J2ni, with ni = Hi{{x{}), (7) 

I 

where the notation in 0 indicates that TLi is to be taken 
as Hi plus an external potential due to the other sub¬ 
systems and represented parametrically. For instance, 
a quantum subsystem acted upon by an external electric 
field is seen in the description of ® as having a Coulomb 
potential energy (operator) coupling all the charges of the 
subsystem with all the charges outside of it which are 
sources of this held, while in the approximate descrip¬ 
tion of 0, it is seen as coupled to an external parameter 
E representing the strength of the held. We shall call the 
latter the thermodynamic description. 

The local equilibrium assumption in the thermody¬ 
namic description is the statement that the macroscopic 
state of each part of our system, with a number of par¬ 
ticles operator Mi, a temperature Ti and a chemical po¬ 
tential p,i, is an equilibrium state. The local equilibrium 


density matrix of the total system is the factorized Gibbs 
state (subsystems macroscopically uncorrelated) 

e' = 0 = (^exp[-(77/ - niMi - fli)/Ti], (8) 

i i 

with Q; = Q,i{Ti,pLi, {;r(^}) the thermodynamic potential 
of subsystem I, introduced so as to normalize the density 
matrix, that is, Tr exp[—{'Hi — fiiMi)/Ti] = exp(-Q;/r;). 
Note that, since the degrees of freedom of different sub¬ 
systems are uncoupled in the thermodynamic description, 
all operators Hi and Mm form a mutually commuting set 
and then dehne a natural basis of common eigenstates 
that we represent as {|a)}. 

The appearance of this natural set defines a family of 
observables acting on the system Hilbert space that, like 
Hi and Mm, we call thermodynamic; these observables 
are diagonal in the basis {|a)}. According to this, 77 is a 
thermodynamic observable, and we denote the set of all 
these operators as 

T={g-.[g,n] = Q} ( 9 ) 

Clearly all constant operators as well as all time-averaged 
observables22iM belong to this family. With V denot¬ 
ing the projection operator to the subspace spanned by 
{|q:)(q;|}, we can then split an arbitrary observable G in 
the convenient form 

G = VG + MG = g + G~, (10) 

where g = HG is the thermodynamically measurable 
part (or thermodynamic part) of G; the complementary 
part being GG = MG = G — g. 

The thermodynamic observables must have the char¬ 
acteristic of being slowly-varying quantitiesi^^— This is 
quantified in our theory by introducing a geometric mea¬ 
sure, A, of how approximate is the thermodynamic de¬ 
scription. For this, let us introduce for an arbitrary ob¬ 
servable B the hermitian operator 

CAB = —i[A,B\ with AgT, (11) 

and consider the simple geometry induced by the Hilbert- 
Schmidt norm ||G^b|| = (Tr G^^^G^s)^/^. It is trivially 
seen that 

l|C'«oll = rGo||=0, (12) 

and by using the Jacobi identity for commutators, to¬ 
gether with [g,'H] = 0, coming from ([HI), we can write 

\\Cnd = WCgcnnl ( 13 ) 

where we identify g = Ggn when the parameters repre¬ 
senting external constraints are fixed in time. Therefore, 
if G-uh tend to the null operator in the norm, i.e. if we 
have A 0 with 

= £o = Eq ^ “ ec,,f\{a\H\a')\'^, 

olol ' 


(14) 
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where H|a) = ea\oL) and eg is the smallest characteristic 
energy in the system (making A dimensionless); then we 
conclude, by using ([ 12 ]), (IT^ as well as the continuity 
of the norm, that the quality of slow variation can be 
expressed as 

11^11 = ll[^,i?]ll=0(A). (15) 

The condition A 0 is physically realized when the 
thermodynamic limit is taken for all the subsystems com¬ 
prising the total system, since in this limit (ff) and (H) 
tend to be indistinguishable for arbitrary states. 

We now give the steps that constitute our general 
method in the next section. Given the density matrix 
pt of the total system, we define the entropy operator as 
the negative of its logarithm, S't = — In pt and, from this 
and the aforementioned discussion, the thermodynamic 
entropy operator as St = 'DSt- Since in the local equilib¬ 
rium approximation the density matrix pt = ff is already 
diagonal in the basis {let)}, we have in this case 

S'" =-Ing'= - piMi-{li), (16) 

' i 


Note that for an arbitrary observable G, the iden¬ 
tity (dG/dt) = d{G)/dt holds whenever the density ma¬ 
trix used to calculate the expectation value satisfies the 
Liouville-von Neumann equation, as is easily proved by 
changing to the Heisenberg picture within the expecta¬ 
tion value operation, where dG/dt = dG/dt — i[G,H], 
with G depending explicitly on time in the Schrbdinger 
picture via the external parameters, and using the known 
identity TvA[B,C\ = TrG[A,Hj. Therefore, as long 
as the local equilibrium density matrix p’’ satisfies the 
Liouville-von Neumann equation, we can commute the 
operation (dG) = d{G) and write (ITOl) as the usual form 
of the first law of thermodynamics. 

The equivalence of (fT^ with the usual form of the first 
law of thermodynamics 

Ti dt{si) = dt{Hi) - PI dt{M) + Y,Fi dtxi, ( 20 ) 

A 

then requires that p^ satisfies the Liouville-von Neumann 
equation idtp’^ = [iL, p*'], where we use the symbol dt as 
a shorthand notation for d/dt. For this to be the case, it 
is necessary from (fT^ that 


where we have used the commutativity of all 'Hi and A/"m 
to express In p’’ = interested in thermo¬ 

dynamic entropy changes as the main observable, then 
the next step is an expression for {dtSt), which we get by 
first differentiating (ITCl) 


dS^' ='^Y “ dpiM - pidSfi - dVti 


Ti 


(17) 


i/Hi — piNi — H;)]. 


Since H; is a function of pi , T/ and of the external param¬ 
eters, implicit in "Hi, we can differentiate the normal¬ 
ization relation Tr exp[—{'Hi — piJ\fi)/Ti] = exp(—flj/T;) 
after variations in these arguments to get, after noting 
that (Gz) = Trp^'Gz = Trp[Gz for a local operator Gz 
acting on the /-th subsystem. 


dni = -J2 Fidx\ - {Afi)dpi - -^{{Hi) - pi{M) - ni), 

(18) 

with FI = —{d'Hi/dx\) being the average force exerted 
by subsystem I on its surroundings to get the displace¬ 
ments dx’’^. Taking expectation value of (13 and substi¬ 
tuting (fT^ we conclude that the average rate of change of 
the total thermodynamic entropy is in this case additive, 
{dS’^) = '^lidS/), with 


Ti {dSD = {dUi) - piidAfi) + Hdxi, (19) 

A 

We have arrived in this way to the first law of thermo¬ 
dynamics, through a line of reasoning originally due to 
Gibbs;^ generalized here to the quantum case. 


\\dtp^ = 0(A), (21) 

since is expressed in terms of thermodynamic observ¬ 
ables. When the thermodynamic limit is taken for each 
subsystem, we have A —A 0, and then the parameters 
should vary with time so slowly that the state p*" can be 
interpreted as “moving” in a locus of equilibrium states, 
so that 11 ftp''11 —A 0 in (|2T|l . These are precisely the qua¬ 
sistatic (or reversible) transformations for which the first 
law involving thermodynamic entropy changes applies, 
hence the superindex “r” standing for reversible, and the 
systematic omission of the time subindex in the variables. 
Note that in this case, the quantity Nc [ 7 ] in (P) vanishes 
for any 7 = {x()^(t), V A, I and t G ts]}- 

A nonzero entropy production appears instead when 
the subsystems are macroscopic at the atomic scale, but 
compared to the size of the total system, they are small 
volume elements, Svi. In this case, an entropy balance 
equation may be obtained from (| 20 l) . by using the rela¬ 
tions 

d,{m) = -Y,Jir+Y, P 2 ) 

ft (M) = - E + E (23) 

which state that the average macroscopic energy and 
number of particles of a given subsystem can only change 
by transport to other subsystems, defining the corre¬ 
sponding currents j!f//' and J///^ in terms of quantities 
proportional to the particle velocities, with an appropri¬ 
ate microscopic account for the heat currents, plus terms 
allowing the technical possiblity of a creation or destruc¬ 
tion of particles induced by the variation of the external 




constrainsts. Substituting these in (l20l) we get 


IV. MASTER EQUATION FOR THE 
THERMODYNAMIC ENTROPY OPERATOR 
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dtiSJ) = 


-T 

Ti V 




dtx{ 


^ m 


(24) 


the first term in the first equality being the entropy pro¬ 
duction term, , and the second one the entropy trans¬ 
port term, . Results consistent with the classical 
theory are obtained when particle creation or destruc¬ 
tion is not observed macroscopically, in which case ((22ll 
and (12311 are just the usual conservation laws (continuity 
equations) and the entropy production in the subsystem 
reduces to the well-known sum of products of thermody¬ 
namic forces times the rate of change of their conjugate 
external parameters 

(classical). (25) 

‘ A 


We generalize the thermodynamic description to in¬ 
clude subsystems which are not distinguished by spatial 
boundaries and which are not necessarily macroscopic at 
the atomic scale. The key point to borrow from ther¬ 
modynamics is the existence of the thermodynamic basis 
{|q:)} and the interpretation of thermodynamic observ¬ 
ables as those which are diagonal in this basis. That is, 
we consider an isolated quantum system (containing pos¬ 
sible reservoirs) which has a Hamiltonian H representing 
the energy of uncoupled subsystems, as before, and study 
the dynamics when the perturbation, V, mixing the de¬ 
grees of freedom of the different subsystems, or a set of 
them, is turned on. 

The Hamiltonian of the total system is then given by 
H = T-L + and the situations of interest include phe¬ 
nomena such as quantum quencheSj^^— or the response 
to applied fields After preparation of the system 

in an initial statistical state of the form 


The presentation given here can be straightforwardly 
generalized by considering local equilibrium Gibbs en¬ 
sembles more general than ([8|), that is, by augmenting the 
thermodynamic entropy operator (I16|) with terms pro¬ 
portional to the components of the macroscopic linear 
and angular momentum operators of each subsystem^^ 
with (|2^ and (l23t expanded to include the conservation 
laws of their respective expectation values. 

Note that we have kept the superindex “r” (although 
not strictly with its original connotation) in because, 
even though the thermodynamic limit is not taken for 
each susbsytem, which would make A —^ 0 and the pro¬ 
cesses necessarily quasistatic, the fact that the volume 
elements, Svi, are macroscopic at the atomic scale still 
implies that A is very small and hence, from (1211) . that 
the variations dtp'^ should correspondingly be very small 
in the norm. As mentioned in section m we then see 
why the classical theory works well for slow processes, 
i.e. those for which the time to get relaxation to equilib¬ 
rium within each volume element is much shorter than 
the time to get equilibrium among them. 

The discussion in this section elucidates the problems 
with the local equilibrium assumption and previous the¬ 
ories of entropy production, which rely on expressions of 
the type (1^01) together with conservation laws, like (1^^ 
and (|23)) . as in the classical theory. As we have made 
explicit, developing a theory of entropy production from 
(|20l) inherently assumes that the correlations among the 
subsystems of a large isolated system are negligible for 
all times, and using (1221) in this theory takes for granted 
that an appropriate mechanical description of the micro¬ 
scopies of heat currents have been univocally achieved. 

We now propose a way to derive an entropy balance 
equation for the subsystems of a general isolated sys¬ 
tem from first principles, starting from the Liouville-von 
Neumann equation for the density matrix of the isolated 
system, which does not rely on the above assumptions. 


po=exp{-So), (26) 

with Sq an arbitrary (in general unbounded) hermitian 
operator with ^ 0, the nonequilibrium state is 

described by the evolved density matrix pt , and we define 
the entropy operator St by 

Pt = exp(-5t), or St = -Inpt, (27) 

which can always be written since the density matrix is 
positive-definite. This exponential representation of the 
density matrix is not new; it is a generalized forn4^4£ 
of the nonequilibrium statistical operator introduced by 
Zubarev 4^4^ and obtained for the case of steady states 
by Hershheld^i. 

As discussed in the previous section, our new thermo¬ 
dynamic entropy operator, St = —'Dlnpt, is obtained 
from St by projecting to the space of operators diag¬ 
onal in the basis {|a)} of eigenstates of Ti. We now 
establish the second law of thermodynamics for nonequi¬ 
librium transformations of the total system. For this, 
we consider for simplicity the specific situation of initial 
states diagonal in the thermodynamic basis, e.g. those 
of local equilibrium form as in ([ 5 ]), for which Sq = 0 
or So = Sq. These initial states are usually assumed in 
practice^i^iiZl, e.g. in transport problems. 

Let us denote the diagonal (or thermodynamic) part 
of the density matrix of the system as 

Qt = F>pt. (28) 

The occupation probability of the state ja) is obtained 
by taking matrix elements Pa-t = (cK|0t|cK)- The proof 
now follows in steps by first using a corollary to Klein’s 
inequality?^ which states that for any concave function 
/(x) we have 


Tt f{gt)>Tr fipt). 


(29) 
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By choosing the concave function f{x) = —a;ln(x), we 
easily get 

- Tr Qt In(et) > -Tr pt Inpt or Sd-t > SyN-t, (30) 
where we have denoted Sd-t = — 

diagonal entrop y^^i^^d^d^ and Sy^-^ is the well-known 
von Neumann entropy. Using the time-invariance of 
SyN;t under the unitary evolution of the isolated system 
together with the fact that the initial state is diagonal, 
so that Sd -,0 = SyN-,Q, then (l30l) implie a^^i^^ 

Sd-,t > Sd-p. (31) 

We use this result and the Husimi-Mori lemma,— 
which states that for any convex function g{x) and state 
!'(/;) we have 


('0l5(Pt)l^) > 5((V'|Pt|^/’)), (32) 

to show that, if we choose the convex function g{x) = 
— ln(x) so that —(a| Inptlo;) > — lnPQ,;t, the thermody¬ 
namic entropy satisfies 

St = {St) = Pa-d{a\ln-pt\a) > Sd-,t > Sq, (33) 


where <So = Sd-p = SyN-,o by the assumption of the initial 
diagonal state. For our isolated system for which there 
is no entropy to be transported outside of its boundaries, 
this proves that St satisfies the second law of thermody¬ 
namics. 

Note that, by splitting pt = Qt + p't und using the 
convenient resolvent representation of the logarithm of 
an operator sum 2 S 


ln(^ + B) = 


dx 


1 


1 


x + A + B 

we can expand the thermodynamic entropy as 


(34) 


= Sd-d 

+ E 

a,/3(#Q;) 


1 Pa-,t . Pf3-t 

XPp-d - Po.-,t) ~ (P/3;t - F’«;t)2 l^t_ 

x|(a|ri/ 3 )P + 0 ((n)^), 

(35) 


with St — Sd-d > 0 due to (1551) : therefore the thermody¬ 
namic entropy, unlike the diagonal entropy, is able to cap¬ 
ture entropy increasing processes due to quantum corre¬ 
lations or entanglement among the different subsystems, 
that are encapsulated in the off-diagonal elements of the 
density matrix. When these quantum correlations are 
negligible which, as discussed in section ITTll is the case 
when each subsystem is macroscopic, the diagonal en¬ 
tropy becomes the thermodynamic entropy according to 
(1351) . and due to the quasistatic (or slow) nature of the 
global transformations involved in this case, the ther¬ 
modynamic basis may be referred to as the adiabatic 
basis— 


The thermodynamic entropy, unlike the diagonal and 
von-Neumann entropies, satisfies the third law of thermo¬ 
dynamics in a transparent way. The third law states that 
the thermodynamic entropy at zero temperature must be 
zero. The standard argument is that at zero tempera¬ 
ture any physical state is pure. For an arbitrary pure 
state l"^), there is always an orthonormal basis of Hilbert 
space which has this state as one of its elements (con¬ 
struct it via the Gram-Schmidt procedure starting from 
I')/')). Denote this basis, {|V'r)}) and order its elements 
such that IV') = IV'i)- We take this basis as the reference 
for “diagonal”. With this we then have for the diagonal 
and von-Neumann entropies 

Sdi'ip) = SyN^ = -y^Pr-lnPr, 

E (36) 

= -1 • ln(l) - ^0 • ln( 0 ), 

r^l 

where Py is the probability that the system be found in 
state IV'r)- Eq. (1551) is usually understood to be zero;^ 
although it is clearly an undetermined quantity since, 
taken at face value, —0 • ln( 0 ) = 0 • oo. 

The thermodynamic entropy of pure states is well- 
defined and readily vanishes. In order to show this, we 
denote the density matrices (projectors) pr = |V’r)(V'r|, 
with 'Yhr Pr = 1- We can then write 

Inpi = ln(i - (37) 

Using this, we can compute the thermodynamic entropy 
of the state |V’) as 

S{'^) = -^(V’rlpiDln(pi)IV'r) = -(V’lln(pi)l'!/')- (38) 


This clearly vanishes exactly since jV') = IV'i) is orthog¬ 
onal to all \'4’r^i) involved in the last equality of (1571) . 
This establishes the third law of thermodynamics. 

We are after an entropy balance equation for the sub¬ 
systems, so we need an equation of motion for St and a 
procedure to get from this one for each subsystems, as in 
the previous section. This can be obtained by first not¬ 
ing that the usual unitary evolution of the density matrix 
implies that St also satisfies the Liouville-von Neumann 
equation's satisfied by pt- We have 

idtSt = [H,St]=LSt. (39) 


This allows us to follow exactly the same procedure origi¬ 
nally used with the density matrix— dS to derive an equa¬ 
tion of motion for its diagonal part, Qt^ the so-called 
Nakajima-Zwanzig generalized master equation. That is, 
we split the entropy operator into a diagonal and non¬ 
diagonal part, with respect to the eigenbasis of as 
St = St + S ^, and obtain an equation of motion for the 
diagonal part using Zwanzig’s integralSS 


idtSt = VLSt + VL e 


— itAfL 


Sq i 


drKySt-r, (40) 
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where the memory kernel is defined as 21 

Kr = VLe-^^^^NL. (41) 

Now, it is easy to verify that VLT) = 0 for any 
Hamiltonianj^i therefore the first term in ((40l) vanishes 
and, with our initial diagonal states implying Sq = 0 , 
we are left with the integro-differential equation 

dtSt = - f drKrSt-r- (42) 

Jo 

Although an exact solution for (H^ . as well as for the 
similar equation satisfied by gt , can easily be found by a 
Laplace transformation followed by an inversion 

I l-c+ioo St 

"«'C>0, (43) 


Moreover, Pa = (Q(|pj°^|a) is the occupation probabil¬ 
ity of the state jcc) in its lowest-order approximation^ 
which also satisfies the Born-Markov limit of the gener¬ 
alized master equation, that is, the transport (or Pauli) 
equation 

dtPc. = Y^^Pc'Wa'a - PcWaa')- (49) 

ol' 

The right hand side of (l47l) can be rearranged to yield 
the quantum version, in the Born-Markov limit, of the 
entropy balance equation. We find 

(ft5t)=n-$, (50) 

where the average rate of entropy produced in the system 
is interpreted as 2 i“— 

n = i E (Po^Waa' - Pc'Wa'a) hr , (51) 

^ ^ot' ct 

OL.Ot' 


where Kg is the Laplace transform of the memory kernel, 
obtained from (HTl) as 

we restrict here, for the sake of a clear presentation and 
for comparison with the classical results, to the Born- 
Markov approximation for weakly coupled subsystems, 
leaving a more general discussion for another publication. 
This approximation, which is justified in the limit of very 
weak coupling potentials, P, and very long times (Van 
Hove limi t^^d^ ) amounts to neglecting memory effects in 
(14211 . In practice, this works for times after any transient 
effect or prethermalization plateaii^Sr— of the isolated 
system has passed. We then have in this limit 


dtSt = - lim KgSt, (45) 

S—>-0+ 

where Kg and 5t, after being expanded in powers of V, 
are truncated up to the lowest orders, for which the well- 
known identity for the resolvent operator expansion 

{A +B)-^ = A-^ - A-^B{A +B)-\ (46) 

is very useful. Taking expectation value of (H5|) . and 
noting that for a diagonal operator Q we have {Q)t = 
Trpt^ = Tr QtG, the average rate of change of the ther¬ 
modynamic entropy in the Born-Markov limit is then 

(dtSt) = Y, PcWac' In (47) 

, ^CX.' 

OtOL 

with the transition rates Waa' = 27rS(ea — £a')|4aa'P, 
calculated in the lowest order in the coupling potential 
using Fermi’s golden rule. Here, we have derived the 
transition rates from (l44l) and (l45|) . by using the repre¬ 
sentation of the delta function^^ 

lim Re—— =7rS(uj). (48) 

s-)' 0 + s + ioj 


and the average entropy flux to the surroundings as 
$ = i E (^c.Waa' - Pc.'Wa'a) In 


Wa'r 


(52) 


Of course, the latter must be zero for an isolated sys¬ 
tem since a global entropy current finds nowhere to go in 
this case. The vanishing of this quantity is clearly seen 
from the symmetry of the transition rates Waa' under 
the interchange of indices, resulting from the hermiticity 
of the perturbation V. A nonvanishing entropy current 
is obtained however when we consider the local entropy 
production in a subsystem of a larger system, as in the 
electrical conduction problem of section ED 

Note that H is a sum of terms of the form (x—y) \n{x/y) 
so is always non-negative. It vanishes for reversible trans¬ 
formations (local equilibrium) or in equilibrium due to 
detailed balance, Waa' = Pf^' Wa'ai this being a sta¬ 
tistical statement of the second law of thermodynamics in 
the Clausius form. The outlined method is the one that 
we shall follow in section ED for the electrical conduction 
problem to derive an entropy balance equation for the 
electronic subsystem in the Born-Markov limit, based on 
the transport equation for the total electrons -I- phonons 
+ field system, without any need to invoke expressions 
like (OUl) together with extra conservation laws. 

One of the advantages of our approach, besides be¬ 
ing grounded on fundamental facts regarding the nature 
of thermodynamic observables is that, as opposed to ac¬ 
tively studied relative-entropy formulations^^— of quan¬ 
tum entropy production, it can be generalized to initial 
states with correlations among the subsystems, i.e. not 
of the local equilibrium form. This is very important 
since the neglection of correlations in the state of an iso¬ 
lated system is inconsistent with the specification of its 
energy— We have safely ignored this fact in our present 
discussion because the consideration of a nonvanishing 
second term in (HOI) , due to Sq 7 ^ 0 , only adds the term 


1 

27ri 


/*C+200 


ds- 


Kg 


-VL- 




(53) 









to the solution (l43l) . However, it is easily seen that ex¬ 
pressions containing Sq contribute higher order terms 
in the weak coupling expansion embodied in the Born- 
Markov limit and then are negligible; the same happens^ 
for the contributions coming from Pq in the Born-Markov 
limit of the generalized master equation for pt- There¬ 
fore, our formalism has room to study memory effects and 
strong correlations in the initial state by only straightfor¬ 
ward modifications. These memory effects are the ones 
responsible for heat transport depending on the path 
of thermodynamic states in phenomenological thermo¬ 
dynamics. 


V. RELATION WITH CLASSICAL 
STOCHASTIC THERMODYNAMICS 


We now show that our result, (03), is consistent with 
the result for the average rate of change of the thermody¬ 
namic entropy obtained in Onsager’s classical theory. We 
consider an isolated macroscopic system which has been 
left alone for a very long time (aged system). The classi¬ 
cal thermodynamic state is described by a set of exten¬ 
sive variables, such as energy, mass, electric charge, etc., 
which randomly fluctuate about their equilibrium values 
and whose values define the classical state of the system. 
This state is represented by the symbol at (shifted to 
vanish in equilibrium), whose successive values in time 
describe a stationary stochastic process. 

It can be shown that, if the fluctuations follow a Gaus¬ 
sian process, which can be argued to be the case if the 
extensive variables are algebraic sums of very many in¬ 
dependent (weakly coupled) “microscopic” quantities so 
that the central limit theorem can be invoked and, if in 
addition the process is Markovian, then the joint proba¬ 
bility distribution, 

B(a', At, a") = Pa'Pa'a" (At), (54) 


for observing the values = a! and aj// = a" at the 
respective times separated by an interval At = t" — t', 
with Pa'a" the corresponding conditional probability 
to make a transition between these states, is given by 
Onsager’s principle^^— which we write as^^ 


2 In f2(a', At, a") = Sa' + Sa" 



where the path of integration is the trajectory, a^-, which 
makes the integral a minimum, subject to the conditions 
at' = a' and at" = a". Clearly, if we take the limit 
At —>• 0, the integral tends to iSa'At, where Sa' is the 
entropy production rate in the state a', whose entropy 
is related to the probability distribution. Pa', by Boltz¬ 
mann’s principle. Subtracting the time-reversed expres¬ 
sion of Onsager’s principle from (I55|) we get, in the limit 
At —^ 0, the alternative form 


, 0(a',At,a") ^ ^ 

"" 0(a", - At, a') ~ 2 


(56) 


We now average (IMl) over the joint distribution (I54L 
which is expanded up to linear order in At by writing 
the transition probabilities to go from a' to a” after a 
time At as 

Pa'a"{^t) = 5a'a" + Wa'a"At = Pa''a'{-^t), (57) 

the last equality being the statement of Onsager’s mi¬ 
croscopic reversibility^iSi and leading to the symmetry 
of the transition rates, Wa'a", under the interchange of 
indices. This symmetry allows to write the averaged left- 
hand side of (HI]) as At X^a'a" Pa'Wa'a" ^a{Pa'/Pa") and 
the right-hand side as At Pa'Sa' Therefore, by rec¬ 
ognizing the latter sum as (5), we get the expression 

{S)=Y,Pa'Wa'a"\a^, (58) 

/ // ^ a" 

a'a" 

which gives the desired link with our theory, by compar¬ 
ing with (I47p . We remark that (1551) is of the same fornt^ 

In 

PAt(-cr) 

as Gallavotti and Cohen fluctuation theoremjS^ if we read 
{l/2){Sa' + ^a"), as a realization of the random number 
a = (\/2){Sa^, +>5a,j„), representing the average entropy 
production in going from af to at" during a time interval 
At along the stochastic trajectory of states; and translate 
the joint probability, B(a', ±At, a"), to have the state re¬ 
alizations at' = a' and at" = a", in a forward (+At) or 
backward (—At) evolution, to the corresponding proba¬ 
bilities PAt(±o') to have the realization, (l/2)(5o'-t-5o"), 
of a or its time-reversed value. 


VI. ENTROPY PRODUCTION IN 
ELECTRICAL CONDUCTION 


We next apply the formalism to a model of indepen¬ 
dent electrons coupled to phonons in the presence of an 
electric field. We are interested in the average rate of en¬ 
tropy produced in the electronic system and transported 
to the phonons in the steady state. The picture is then 
that of a large system divided into three subsystems, the 
electrons, the phonons, and the sources of the held. In 
the thermodynamic description we parametrize, as usual, 
the coupling to the latter by introducing Et and forget¬ 
ting about the structure of this subsystem. 

The Hamiltonian of the total system is then 

H = Hei + Hph + Hei-ph + Hp-t, (60) 

where He\ = cjjcfc is the kinetic energy opera¬ 

tor for the electrons, which are assumed to be free ex¬ 
cept for their interaction with the held and the phonons, 
the energy operator of the phonon subsystem is i7ph = 
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'Yhq^q the electron-phonoii interaction is bilin¬ 

ear in electron operators and linear in phonon operators 

He\-ph = 'y ] ^k'k (d? + ( 61 ) 

qkk' 

with representing the strength of the coupling. The 
generalization to multiple electronic bands and multiple 
phonon branches is straightforward and does not change 
the results. Finally, HF-,t represents the effects of the 
applied electric field, Et, and can be written in first- 
quantized notation as HF-,t = -eEt- where Xe 

is the displacement of electron e from some arbitrarily 
chosen reference position. 

Up to time t = 0 we have a collection of electrons in 
local equilibrium with the lattice vibrations of a metal 
at a temperature T, and no applied electric electric field, 
i.e. Eo = 0. The initial state is then of the form 

Po = Z~^ exp[-(?io - /r ATeO/T], (62) 

where Z = ZgiZpi^ is the grand partition function, A/^ei is 
the operator for the total number of electrons, and TLo is 
the Hamiltonian of the uncoupled subsystems 

?io=-ffel + 4h, (63) 

whose eigenstates, constituing the thermodynamic basis, 
are 

|a) = |nin 2 ■ • • Ufc • • •)|iViiV 2 ■ ■ ■ Ng ■ ■ ■) = \n; N), (64) 

which represent the number of electrons, {uk}, and 
phonons, {A(j}, in each single-particle state. 

The electric field is turned on at time t = O’*' to a con¬ 
stant value, i.e. Et = E ioit > 0, and the subsystems are 
subsequently coupled. In the notation of section HVl we 
then have in the generalized thermodynamic description 

H = U = I?el-ph. (65) 

Note that V is the coupling which fully mixes the de¬ 
grees of freedom of the different subsystems (like the Jlim 
in section EB, which need not be separated by spatial 
boundaries. We now explain with some detail how the 
perturbation scheme developed in section IIVI applies to 
the present case. However, we only need to concentrate 
on how the transport equation is obtained in the Born- 
Markov limit, since this suffices to get the average rate 
of entropy production. 

The idea is then to first derive the transport equa¬ 
tion for the total system from the Liouville-von Neumann 
equation; we do it much in the same spirit as Kohn and 
Luttinger— did for elastic electronic scattering and gen¬ 
eralized by Argyres^l to inelastic scattering. Having this 
transport equation, the average rate of change of the to¬ 
tal thermodynamic entropy in the Born-Markov limit is 

=-^(atP„)lnP„, (66) 


as can easily be verified by using (l49l) in (|Tf)) . By proceed¬ 
ing with the transport or quantum Boltzmann equation 
for the electronic subsystem, we obtain a simple expres¬ 
sion for the electronic entropy production. 

For the purpose of the present discussion, it suffices to 
work with the Liouville-von Neumann equation to first 
order in the electric field. That is, with pt = po + Pi;t 
and linear in the electric field, we write 

idth-,t = [Po + V, pi-t] + [Hf, Po], (67) 

where pi;o = 0. The Laplace transform of this equation, 
with pi;s = e“®‘pi; 4 , reads 

ispi-s = (Uo -I- Lv)pi-,s + s~^LfPo- (68) 

With pi-s = Ppi;s and p^.^ = Afpi-s, we separate this 
equation into a diagonal and a non-diagonal part, ob¬ 
taining, respectively, the coupled algebraic equations 

ispi-s = E>Lv pZ,s + s~^T^LfPo, (69) 

[is -l- Af{C q + Ly)] Pl.g = AfLy pi-^s s ^A/LfPo- (70) 

Solving for pi.^ in (1701) and substituting the result in (1001) 
we get a decoupled equation for pi-s, which in the lowest 
Born approximation for the electron-phonon scattering 
reads 

isQi-s = VLy- - MLy Qi-s + s~^VLfPq. (71) 

is + Al Lq 

From this, the transport equation for the total system 
easily arises in the Born-Markov limit by taking the 
Laplace inverse and neglecting memory terms. In terms 
of the occupation probabilities Pa = (ajptja) we get 

dtPa = ^{LfMc + Y.^P<.'Wa'cc - Pc^Wao.') (72) 

a' 


with the transition rates induced by the electron-phonon 
coupling Waa' = 27ri5(£:ct—Ect')IIPei-phIciOP- We have 
then derived the transport equation for the total system, 
in terms of which the average rate of change of the total 
thermodynamic entropy can be calculated, in the Born- 
Markov limit, using (155)) . 

To proceed with the calculation of the entropy produc¬ 
tion of the electronic subsystem, we note that 

p _ pel pph el-ph ('7Q\ 

-‘a ~ -‘AT XnN ’ V ' 

where P®* is the probability that the electrons are in the 
Fock state |n) regardless of the state of the phonons, P^ 
is the probability that the phonons are in the Fock state 
|iV) regardless of the state of the electrons, and is 

the conditional probability that the total system is in the 
state la) in ()64l) . given that the electron and phonon sub¬ 
systems are in states |n) and \N), respectively, without 
“knowing” about each other. Clearly, XnN^ i® ^ function 
of the electron-phonon coupling strength, and can then 
be expanded in a power series of it 


el-ph 1 , el-ph(l) 

XnN = 1 + XnAT 


el-ph(2) 

■XnN 


(74) 
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In the lowest Born approximation for the electron- 
phonon scattering, the electron and phonon subsystems 
are uncorrelated, i.e. = 1, which is the usual Born- 

Oppenheimer approximation, and then by substituting 
(172)1 and (173)1 into (IMll , the average rate of change of the 
thermodynamic entropy of the total system turns out to 
be additive. For the electronic subsystem we have 

(a*5t)ei = -^(9tP:')lnP,?', (75) 

n 

where the normalization condition P’^ = 1 has been 
used. Here, the transport equation for the electronic sub¬ 
system is obtained from (17211 by summing over N 

dtPn = -Y.^LFMr.N,uN + Y,i.Pn'^n’n " Pn^n^'), 

* N 

(76) 

where we have dehned the phonon-averaged reduced 
transition rates r„„/ as 

Tnn' =J2 Pn''Y. WnN,n'N' ■ (77) 

N N' 

We can still go further and use the assumed statistical 
independence of the electrons to factorize their probabil¬ 
ity distribution into the probabilities of the one-electron 
states 

Pn ~ PniPn2 ' ' ' Prik ' ' ' ^ 

where Pn^ is the probability that the one-electron state 
with quantum number k has occupation Uk = 0,1. Sub¬ 
stituting this in (17511 we obtain an additive contribution 
to the average rate of change of the thermodynamic en¬ 
tropy of the electronic subsystem 

{dtSt)el = -^{dtPn^)\np„^ = 

k,nk k 

(79) 

where in the last equality we identify the nonequilibrium 
one-electron distribution as fk = '^kPu^ = Pnk=i 
and use = 1 to express Pnk=o = ^ - fk- The 

transport equation for fk is obtained by multiplying (l76ll 
by Uk and summing over all n. To this end, note that 

Fnn^ — ^ ^ '^kk' ^fc(l ^fc') 

k,k' {k^k') (80) 

X |(- ■ - nfe/ - 1 • • - rife -b 1 ■ • • |n')p, 

which is obtain by using (16111 explicitly, where the one- 
electron transition rate from state k to state k' is 


number q. We assume that the phonon subsystem can be 
kept in equilibrium at temperature T (hence the depen¬ 
dence of N on LOq only), no matter the nonequilibrium 
state of the electrons, as is the case for a good enough 
heat reservoir. 

That the phonons can be considered as a heat reservoir 
in the Born-Oppenheimer approximation can be seen by 
looking at the transport equation for the phonon subsys¬ 
tem, obtained from (17211 by summing over n 

dtPft = - PftQNN'), (82) 

N' 

where we have defined the electron-averaged reduced 
transition rates as 

©WAT' = ^2 Pn ^2 ^nN,n'N'- (83) 

n n' 

Here we observe the important fact that the contribution 
from the first term of (17^ vanishes due to the null value 
of the trace of the commutator [Hf-, Po] in the subspace 
of electrons. This allows the existence of a steady state 
solution of (15711 for which detailed balance holds, which 
is then an equilibrium solution. In any case, the assump¬ 
tion that the phonons are in equilibrium is not necessary 
for the following derivation of the electronic entropy pro¬ 
duction, as N{u}q) in ()5T)l can be replaced by the more 
complicated average obtained by using the nonequilib¬ 
rium solution of (182)1 . not investigated here. 

The transport equation for the one-electron distribu¬ 
tion is then found to be, from (175)1 

dtfk = - ^ nfe(Li7’po)„Ar,„iv+^ P®'r„„/). 

nN nn' 

(84) 

This is just the quantum Boltzmann equation. To write 
it in the familiar form we first note that, by writing H^i 
in first-quantized form and using the well-known formula 
[Xe,f{pe)] = P^pJiPe) we have LfPo = i{e/T)E ■ vpo, 
where i) = '^^{Pe/'kn) = J^k '^kcl.Ck is the velocity oper- 
ator of all electrons, with Vk = the band velocity. 

Therefore, the first term in (f 8 ^ is 

/fc)drift — V ^ ^ ‘^k{^FPo}nN,nN ~ ~Tfr ’ i^k'^ Po) ■> 

I 1 

nN 

= eE • Vkfk{l - fk)/T = -eE ■ Vfc/°, 

(85) 

where /)? = Tr(hfe/5o) is the equilibrium Fermi-Dirac one- 
electron distribution. The second term in (IMll can be 
written, using (ISOll . as 


Wkk' N{ujq)S {ek' - Cfc - uiq) 


-|- [l -|- iV(Wq,)] 6 {ek' — Cfc + Wq) 


( 81 ) 


with N{ujq) = J2n P^i^\^q^q\^) average num¬ 
ber of phonons in the single-particle state with quantum 


{dtfkUn = 

nn' 

= '^P^nk y] Wk'k"{^ - nk')nk" 

n k',k" 

- ^ Punk y] Wk'k"nk'{l - Uk"). 

n k',k" 


(86) 
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Therefore, by noting that nk{^ — Uk) =0 and using (l78l) 
we see that the terms which do not cancel in the above 
sums are 

k)co\\ ~ 'y '^k'k (1 fk') fk'^kk' (1 fk')]- 

k' 

(87) 

We have thus arrived to the familiar form of the quan¬ 
tum Boltzmann equation by substituting (IM|) and (|^ 
into (15^ . With this, we can rewrite the average rate of 
change of the thermodynamic entropy of the electronic 
subsystem, from (17^ . as 

{dtSt)el = Hel - $el, (88) 


which is the entropy balance equation for the electronic 
subsystem, with the average electronic entropy produc¬ 
tion rate 


IIcl 


^ (1 “ (1 “ •Cfe)] 


. fk' Wk'k (1 — fk) 
"" '' fkWkk'{l- f'k)^ 


(89) 


which, similar to (EB, is a sum of terms of the form 
{x — y) hY{x/y) and then satisfies the second law of ther¬ 
modynamics; and the entropy flux from the electrons to 
the phonons is 




kk' 


^{dtfk)dnft In 


fk 


I-fk 


Wk'k 

Wkk' 

(90) 


In the steady state the left-hand side of (1551) is exactly 
zero and then all the entropy produced in the electronic 
system is transported to the phonons. We now want 
to show that this steady-state entropy flux toward the 
lattice vibrations gives the expression of the well-known 
Joule heating. 

We need a solution, /^ = /)) -|- Sft, of the quantum 
Boltzmann equation which, to linear order in the electric 
field strength, we write formally as 


k' 


eE-Vk' ,0 


T 


fk'{l-fk') 


(91) 


where the linearized collision operator W, has matrix el¬ 
ements 

Wkk' = fkWkk' + Wk'k (1 - fk) - — ■ (92) 

'^k 

The quasiparticle relaxation time Tk is given by 

and becomes equal to the momentum relaxation time if 
the transition rates Wkk' are independent of the angle 
between k and k'. 


We now expand (15^ . Because both the logarithm and 
the prefactor vanish in equilibrium the leading contribu¬ 
tion is 0(5f)^. The term from expanding the logarithm 
is easily seen to be 


5fk' 


5fk 


fk'i^-f°k') f°k{^-ny 

while the term coming from the prefactor is 
Wkk'Sfk' - Wk'kSfk- 

Combining these equations with (1911) yields 

n.. = Z 


T 


T 


(94) 


(95) 


where, in the last equality, we recognize the electric cur¬ 
rent as e{v) = eVkSfk = crE , with a the electric 
conductivity. Thus we see that to leading order in the 
electron-phonon coupling and the electric field, and on 
the assumption that the phonons act as a reservoir, the 
electronic entropy production predicted by our formula 
is exactly the result expected from the Joule heating, 
TIIoi = aE^, implied by the electric field. Therefore, 
as desired, we have arrived at an expression of energy 
dissipation from a first-principle calculation of entropy 
production, not the other way around, as in previous ap¬ 
proaches. 

We remark that the results for the entropy produc¬ 
tion presented here are beyond the linear response the¬ 
ory. This is because, even when starting from the linear 
in the electric field correction to the density matrix, pi-t 
(see Eq. (l67l) l. we derived the leading contribution to the 
electronic entropy production which is quadratic in the 
electric field. This is in contrast to past approache a^^i^^ 
for the calculation of the Joule heating, which requires 
going to the second order in the electric field contribution 
to the density matrix p 2 -,t for the calculation of the rate 
of change of the energy of the electrons. A field-theoretic 
approach^SiSS for the calculation of higher order terms 
in the entropy production, beyond the Born-Markov ap¬ 
proximation will be treated elsewhere. 

It is illustrative to evaluate the result explicitly, assum¬ 
ing e.g. dispersionless optical phonons uiq = wq. With 
l-^fe'fcP = M 5q^k'-k, and assuming a degenerate electron 
system (i.e. T -C ep) we obtain 


Ilel 


Eep—ujQ T E 


{ep/T) sinh(a;o/r). 


(96) 


with the electron density of states. In this case, the 
entropy production becomes large at low temperatures 
due to an increase in the conductivity (phonons not ther¬ 
mally activated and then scarcity of scattering centers), 
and hence in the Joule heating; this is expected when the 
only scattering mechanism is from optical phonons. 

Finally, we would like to point out the connection of the 
result (1551) with the discussion in section m concerning 
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the foundations of the classical theory. With only the 
action of one of the subsystems (the sources of the E- 
field) treated parametrically, with the three spatial com¬ 
ponents of the field E\ playing the role of the external 
parameters to the electronic subsystem, we can define an 
operator F\ = dHp/dEx for the force exerted on the 
electrons upon variation of the field and write 

Tflei = X] FxdtEx = dtJ2 F^Fx - J2idth)Ex 

A A A (97) 

= dtHp + ev ■ E = ei) ■ E, 

where Hp = ^Exx^, and to get the last equality 

we use dtHp = dtH = 0, since the total system is iso¬ 
lated. Taking expectation value of (IWll . the last equality 
is just (EH) and the form of the first equality is reminis¬ 
cent of the classical expression (1251) . 

We then see that, although in the present discussion 
the subsystems are not separated by spatial boundaries 
(the essence of the generalized thermodynamic descrip¬ 
tion) and there is no local equilibrium at all times: the 
phonons remain in equilibrium, as implied by the as¬ 
sumption that they constitute a good heat reservoir, but 
the electrons attain a nonequilibrium steady state; the 
common feature with the discussion in section Hill is the 
complete factorization of the probability distribution of 
the system over the degrees of freedom of the different 
subsystems (uncorrelated subsystems), here manifested 


as the Born-Oppenheimer approximation. An appropri¬ 
ate account of the quantum correlations between subsys¬ 
tems is therefore the key to purely quantum thermody¬ 
namic behavior. 


VII. CONCLUSION 

We have developed a theory for the entropy production 
in quantum many-body systems by introducing an en¬ 
tropy operator and calculating the average rate of change 
of its thermodynamically measurable part. We show that 
the laws of thermodynamics are satisfied exactly within 
our formalism. In the Born-Markov approximation which 
describes the physics of weakly-coupled subsystems of 
an isolated system in the long-time limit, the theory re¬ 
produces the entropy balance equation which is funda¬ 
mental in classical nonequilibrium thermodynamics and 
the Joule heating contribution to the entropy produc¬ 
tion expected in a standard conductor. Applications to 
other systems as well as generalizations beyond the weak- 
coupling limit will be presented elsewhere. 
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